Half-Metallic Graphene Nanoribbons 
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Electrical current can be completely spin polar- 
ized in a class of materials known as half-metals, 
as a result of the coexistence of metallic nature for 
electrons with one spin orientation and insulating 
for electrons with the other. Such asymmetric 
electronic states for the different spins have been 
predicted for some ferromagnetic metals— for ex- 
ample, the Heusler compounds-*^ and were first 
observed in a manganese perovskite^. In view 
of the potential for use of this property in realiz- 
ing spin-based electronics, substantial efforts have 
been made to search for half-metallic materials^'*. 
However, organic materials have hardly been in- 
vestigated in this context even though carbon- 
based nanostructures hold significant promise for 
future electronic device®. Here we predict half- 
metallicity in nanometre-scale graphene ribbons 
by using first-principles calculations. We show 
that this phenomenon is realizable if in-plane ho- 
mogeneous electric fields are applied across the 
zigzag-shaped edges of the graphene nanoribbons, 
and that their magnetic property can be con- 
trolled by the external electric fields. The results 
are not only of scientific interests in the interplay 
between electric fields and electronic spin degree 
of freedom in solids^ but may also open a new 
path to explore spintronics^ at nanometre scale, 
based on graphene*'^'-'^*''-'^-'^ 

When a single graphite layer is terminated by zigzag 
edges on both sides, which we refer here to as a zigzag 
graphene nanoribbon (ZGNR) (Fig. 1), there are pecu- 
har locahzed electronic states at each edgei^i^^. These 
edge states (which are extended along the edge direc- 
tion) decay exponentially into the center of the ribbon, 
with decay rates depending on their momentu m^^i^'^i^^i^^ . 
Such states have been observed in monoatomic step edges 
of graphite by using scanning probe techniqueai^iii. The 
localized edge states form a two-fold degenerate flat band 
at the Fermi energy [Ep], existing in about one third of 
the Brillouin zone away from the zone center-i^JJiiiiii^. 
By invoking band ferromagnetism, it has been suggested 
that an opposite spin orientation across the ribbon be- 
tween ferromagnetically ordered localized edge states on 
each edge in ZGNRs is the ground-state spin configura- 
tion; that is, the total spin is zero ^^'^^i^^ . Because the 
states around Ep are the edge states and linear combina- 
tions of them, the effects of external transverse fields are 
expected to be significant on these states, in constrast 
with those on the extended states2£. 

Our study of the spin resolved electronic structure of 



ZGNRs is based on the ah initio pseudopotential density 
functional method^^ within the local spin density approx- 
imation^^. A periodic saw-tooth- type potential perpen- 
dicular to the direction of the ribbon edge is used to sim- 
ulate the external electric fields (-Ecxt) in a supercell (Fig. 
1). In accordance with previous conventiorki^ii^iiiii^, the 
ZGNRs are classified by the number of zigzag chains (n; 
Fig. 1) forming the width of the ribbon. We will here- 
after refer to an ZGNR with n zigzag chains an n-ZGNR. 
When the spin degree of freedom is neglected, our calcu- 
lation from first-principles also predicts a two-fold degen- 
erate flat band at Ep (Fig. 2a). But the spinless state is 
not the ground state. Moreover, the electronic structures 
of the ZGNRs show marked alterations when spins and 
i?cxt are included. 

Considering first the spin degree of freedom, we find as 
in previous studies that the configuration with opposite 
spin (antiferromagnetic) orientation between ferromag- 
netically ordered edge states at each edge (Fig. 2b) is 
favored as the ground state over the configuration with 
same spin orientation between the two edge o-^^'^^i^^ . (The 
present result of antiferromagnetic spin configuration on 
the honeycomb lattice is consistent with a theorem for 
electrons on bipartite lattice^^.) Our calculations show 
that the magnetic interaction energies are quite large. 
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FIG. 1: Graphene nanoribbon in electric fields. Dia- 
gram of a zigzag graphene nanoribbon (ZGNR) with exter- 
nal transverse electric field -Bcxt. £cxt is applied across the 
ZGNR along the lateral direction [x direction) in an open cir- 
cuit split-gate configuration and is positive towards the right 
side. The ZGNRs are assumed to be infinite along the y direc- 
tion. A small longitudinal source-drain field could be applied 
to generate spin-polarized currents along the {/-direction. Hy- 
drogen atoms on the edges are denoted by circles. The ZGNR 
shown in this figure is an 16-ZGNR (n=16). The width w of 
ZGNRs under study was in the range 1.5 ~ 6.7 nm. 
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For example, the total energy difference between a spin- 
polarized edge and a spin-unpolarized one is 20 meV per 
edge atoms in the case of 8-ZGNR, and the spin config- 
uration is further stabilized by 2.0 meV per edge atom 
due to the antiferromagnetic coupling between the spin- 
polarized edges. Because the interaction between spins 
on opposite edges increases with decreasing width, the to- 
tal energy of an n-ZGNR with antiferromagnetic arrange- 
ment across opposite edges is always lower than that of a 
ferromagnetic arrangement if n < 32. This total energy 
hierarchy is maintained when external electric fields are 
applied. It is known that spontaneous magnetic order- 
ings in one- and two-dimensional spin lattice model are 
difficult to achieve at finite temperature^^. Spin correla- 
tion lengths comparable to nanoscale systems, however, 
is possible in practice^s^i^liS^. Here, we also expect 
that spin orderings are realizable because of the large 
anisotropic exchange interactions between the spins in 
ribbons with split-gate geometry on the substrate. 

We find that the ground state of the ZGNRs, including 
spin degree of freedom, has a bandgap which is inversely 
proportional to the ribbon width. However, the energy 
splitting at fca = TT is ~0.52 eV regardless of width if n > 
8. The states of opposite spin orientation are degenerate 
in all bands (Fig. 2c, left). When spins are included, the 
degeneracy between the occupied and unoccupied edge- 
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FIG. 2: Electronic structures of graphene nanorib- 
bons. In all the figures, the Fermi energy (Ef) is set to zero, 
a. The spin-unpolarized band structure of an 16-ZGNR. b. 
The spatial distribution of the charge difference between a- 
and /3-spin {pa{r) — Ppif)) for the ground state when there 
is no external field. The magnetization per edge atom for 
each spin on each sublattice is 0.43 ps with opposite orien- 
tation, where ps is the Bohr magneton. The graph is the 
density integrated in z direction and the scale bar is in unit 
of 10~^|e|/A^. c, From left to right, the spin-resolved band 
structures of an 16-ZGNR with E^^t = 0.0, 0.05, and 0.1 
V/A, respectively. The red and blue lines denote bands of a- 
spin and /3-spin states, respectively. Inset, the band structure 
with i?ext = 0.1 V/A in the range \E — Ef\ < 50 meV and 
0.77r < ka < n (the horizontal line is Ep)- 



state bands at Ep is now lifted and the edge states near 
Ep have dispersion along the direction of the edge with 
a band width of ^2 eV when extended over the Brillouin 
zone. 

With applied transverse electric fields, we find that the 
valence and conduction edge-state bands associated with 
one spin orientation close their gap, whereas those asso- 
ciated with the other widen theirs (Fig. 2c). So, under 
appropriate field strengths, the ZGNRs are forced into a 
half-metallic state by the applied electric field, resulting 
in insulating behavior for one spin and metallic behavior 
for the other. We shall defer the discussion of spin-orbit 
interactions later. For now, we label the gap-opening 
states as a-spin (shown in red in Figs 2-4) and the gap- 
narrowing states as /3-spin (blue). In a 16-ZGNR, the 
bandgap associated with /3-spin is completely closed by 
an i^cxt = 0.1 V/A, whereas the gap for a-spin electrons 
remains very large at 0.30 eV (Fig. 2c). The energy gap 
for the /3-spin electrons changes to an indirect gap from 
a direct gap as E^xt increased, and is closed indirectly 
(Fig. 2c, inset). After gap closure, an electron channel 
near ka — n and a hole channel near ka — O.TStt appear 
at Ep, all with the same spin direction. 

The half-metallicity of the ZGNRs originates from the 
fact that the applied electric fields induce energy level 
shifts of opposite signs for the spatially separated spin 
ordered edge states. Such separate and opposite energy 
shifts are made possible by the localized nature of the 
edge states around Ep. Because oppositely oriented spin 
states are located at the opposite sides of the ZGNR, 
the effect of i?ext on them is opposite, moving the occu- 
pied and unoccupied /3-spin states closer in energy but 
moving the occupied and unoccupied a-spin states apart 
(Fig. 3). The electrostatic potential is raised on the right 
side and lowered on the left side as £^cxt(> 0) increases. 
Correspondingly, the energies for localized edge states 
on the right side are shifted upward and those on the 
left side downwards, eventually leaving states of only one 
spin orientation at Ep (Fig. 3b). In a 16-ZGNR, the oc- 
cupied a-spin states and unoccupied /3-spin states on the 
left side move downwards in energy by 19 meV and 110 
meV, respectively, and occupied /3-spin and unoccupied 
a-spin states on the right side upwards by 112 meV and 
74 meV, respectively, as i^ext increases to 0.1 V/A (Fig. 
3c). The occupied /3-spin states in the middle of a 16- 
ZGNR are the tails of the localized /3-spin states on the 
right side and the unoccupied /3-spin states are from the 
left side, so that occupied and unoccupied /3-spin states 
in the middle of the ZGNR move oppositely to close the 
gap. The energies of the occupied and unoccupied a-spin 
states in the middle also follow movements of those of the 
corresponding localized states on each side, resulting in 
an increased gap value (Fig. 3c). 

We note that the critical electric field for acheiving 
half-metallicity in ZGNRs decreases as the width in- 
creases because the electrostatic potential difference be- 
tween the two edges is proportional to the system size. 
For a 32-ZGNR whose width is 67.2A, E^^t = 0.045 
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V/A is required to close the band gap for the /3-spin 
electrons (Fig. 4). Because the energy shifts of the edge 
states depends on the total voltage drop between the two 
sides, the variation of the energy gap is expected to ex- 
hibit a universal behavior as function of wEcxt where w 
is the width of the ZGNR. This is seen in the inset in Fig. 
4. From the calculations, the required critical field is es- 
timated to be 3.0 {Y)/w{k). To establish half-metallicity, 
the relevant energy scale is given by the field-induced en- 
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FIG. 3: Origin of half-metallicity. a, Schematic den- 
sity of states diagram of the electronic states of a ZGNR in 
the absence of an applied electric field. Top: the occupied 
and unoccupied localized edge states on left side (L-region as 
defined at the bottom) are a-spin and /3-spin states, respec- 
tively, and and vice versa on right side (R-region) with the 
same energy gap for both sides. Bottom: schematic diagram 
of the spatial spin distribution of the highest occupied valence 
band states without an external electric field, b, Top: with 
a transverse electric field, the electrostatic potential on the 
left side is lowered {eAV < 0) while the one on the right side 
raised {eAV > 0). Correspondingly, the energies of the lo- 
calized state at the left edge are decreased and those of the 
localized states at right edge increased. Bottom: the result- 
ing states at Ep are only /3-spin. c, From left to right, the 
local density of states of a- and /3-spins (ordinate) of a 16- 
ZGNR as a function of energy (abscissa) for atoms in the L, M 
and R regions shown in a, respectively. From top to bottom, 
Eext= 0.0, 0.05 and 0.1 V/A respectively. The local density 
of states in the middle panels are enlarged eightfold for clar- 
ity For Sext =0.1 V/A, the van Hove singularities of /3-spin 
in the M and R region are above the Ep by 5 meV, and all 
states at Ep are of /3-spin. 



ergy shift, and its magnitude is in the order of 100 meV. 
Thus the small magnitude of spin-orbit interaction(4~6 
meV) in carbon atom o^^'^° would not change the half- 
metallic nature of the ZGNRs, but would function in 
determining the spatial direction (normal direction with 
respect to the ribbon plane) of spin up and down in the 
ZGNR&29. 

Because edges are inevitably susceptible to defects, 
we have examined the robustness of the predicted half- 
metallicity to edge defects. Our calculations show that 
the system remains purely of one spin- type at -Ej^ in the 
presence of different type and concentration of defects. 
Results on 8-ZGNRs with three different kinds of de- 
fects (dangling bonds, vacancies, and Stone- Wales defects 
at 6^12% defect concentration per edge) are presented 
in Supplementary Fig. 1, confirming that the predicted 
half-metallicity is indeed robust. 

Another consideration is that, when in the half- 
metallic state the ZGNRs are in a transverse electric field, 
the current-carrying electrons moving from the source 
to the drain in the longitudinal direction would experi- 
ence an effective magnetic field due to spin-orbit inter- 
actions^ii and the spins are expected to rotate. How- 
ever, we find that the resulting extremely weak effective 
magnetic field are paralle to the spatial spin direction (z- 
dircction in Fig. 1) already determined by the instrinsic 
spin-orbit interaction of carbon atoms. Supposed that 
we have /3-spin electrons moving with velocity v = vy in 
E = i?cxti, the effective magnetic field exerted on the 



S- 0.3 

CO. 







1 


' n=8 • • 




' •n=1100 




n=16B ■ ~ 




n=32A A 


\ \\ \ ° 


c 










\ VSl V 2 






w-E,„(V) 




1 1 



0.1 0.2 0.3 

Eext(V/A) 

FIG. 4: Dependence of half-metallicity on system size. 

Ac (red) denotes the direct band gap of a-spin, and 
(blue) the (in)direct gap of /3-spin as function of _Bext for 
8-ZGNR (filled circles), 11-ZGNR (open circles), 16-ZGNR 
(squares), and 32-ZGNR (triangles). The slope variation 
of Aq and A^ indicates gap change from direct to indi- 
rect. The rescaled gap 5a:{w, -Bcxt) = Adw, _Ecxt)/Ao(ui) and 
Sp{w, Ecxt) = A/3(w, i?ext)/Ao(u)) for the various widths col- 
lapse to a single function of W-Eext as shown in the inset, where 
Aa{w, Eext) and A/3(w;,£'ext) are the bandgap for a-spins 
and /3-spins, respectively, of the ZGNR with a width of w in 
Sext(/ 0), and Aoiw) = Aciw,E^^t = 0) = A^(«;,£;ext = 0). 
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/3-spin electrons would be BcS = -^^Ecxtz where h is 
the Planck constant, m is the mass of an electron, e is 
the charge on an electron, and c is the speed of light. At 
a critical electric field of 0.045 V/A for a 32-ZGNR, the 
estimated energy for spin-orbit coupling due to ii'oxt is 
only 1.1 X 10^"* meV. We also find that, as a result of the 
energy gap asymmetry for each spin, there is no spin pre- 
cession even when the direction of -Ecxt is tilted or when 
the spatial spin direction is altered by spin-orbit interac- 
tion arising from the substrate. So,the spatial spin direc- 



tion once determined would not change even if a strong 
transverse electric field is applied. This implies that, if 
we change the direction of i?ext, the spin-polarity of the 
carriers at Ep of the half-metallic ribbon will be reversed 
because the induced potentials at the edges change their 
signs. Hence, under these conditions, the half-metallic 
nature is robust even though a transverse electric fields 
is applied, and spin polarized current should be obtained 
in transport experiment with split-gates. 
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FIG. 5: Supplementary Figl. Robustness of half- 
metallicity in defective graphene nanoribbons. a, A 

ball-and-stick model of an 8-ZGNR with one dangling bond 
(an edge carbon atom without hydrogen passivation) on the 
right edge. One dangling bond per 17 edge atoms are consid- 
ered, which corresponds to 5.9% impurity concentration per 
one edge. In the figure, the atomic structure near the defect 
is displayed. Electric fields (Scxt > 0) are applied from the 
left side to the right side. The Q-spin state is located on the 
left side and the /9-spin state on the right side in the case 
of the ground state without applied electric fields, b. From 
top to bottom panels, the spin resolved total density of states 
(TDOS) are drawn for a defective 8-ZGNR shown in a with 
Ecxt ~ 0.0, 0.1, and 0.2 V/A respectively. At the same critical 
Eact of 0.2 V/A for an ideal 8-ZGNR, the gap for /3-spin state 
is completely closed, c, A ball-and-stick model of an 8-ZGNR 
with one carbon-atom vacancy on the right edge. Defect con- 
centration, ground spin configuration, and direction of electric 
fields are identical to the case of an 8-ZGNR with one dan- 
gfing bond shown in a. d. The TDOS for an 8-ZGNR with 
one vacancy on the right edge (shown in c) with and without 
electric fields. Half-metallic nature persists also in this case, 
e. A ball-and-stick model for an 8-ZGNR with one single ro- 
tated bond (Stone- Wales defect at the right edge). Impurity 
concentration per one edge in this case is 11.8% since two 
edge atoms participate to create a Stone- Wales (SW) defect. 
Electric fields and ground spin configuration follow the same 
convention described in a and c. The total energy of an 8- 
ZGNR with a SW defect is much higher than that of the ideal 
one by 3.05 eV per defect so that appropriate treatments (e.g. 
annealing) on the sample will remove this highly unstable de- 
fect of this kind. Nevertheless, even with such high impurity 
concentration shown in e, the TDOSs displayed in f clearly 
show the robustness of the half-metallicity since the gap for 
the /3-spin states is completely closed while the gap for a-spin 
state is at 0.39 eV with Ecxt = 0.2 V/A. 



